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ABSTRACT 

We present an analysis of a 1 60 ks NuSTAR observation of the nearby bright Seyfert galaxy IC 4329A. The high- 
quality broadband spectrum enables us to separate the effects of distant reflection from the direct coronal continuum, 
and to therefore accurately measure the high-energy cutoff to be E cut = 1 78 +7 i j ) keV. The coronal emission arises 
from accretion disk photons Compton up-scattered by a thermal plasma, with the spectral index and cutoff being due 
to a combination of the finite plasma temperature and optical depth. Applying standard Comptonization models, 
we measure both physical properties independently using the best signal to noise obtained to date in an active 
galactic nucleus over the 3-79 keV band. We derive kT e = 37 47 6 keV with r = 1 . 25 +( , ) j 2 [ ( ( ) ) assuming a slab geometry 
for the plasma, and kT e = 332°, keV with r = 3.4 1 4 q 5 3 8 for a spherical geometry, with both having an equivalent 
goodness-of-fit. 

Key words: accretion, accretion disks - galaxies: active - galaxies: individual (IC 4329A) - galaxies: nuclei - 
galaxies: Seyfert - X-rays: galaxies 
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1. INTRODUCTION 

The primary hard X-ray continuum in Seyfert galaxies arises 

from repeated Compton up-scattering of UV/soft X-ray ac- 

cretion disk photons in a hot, trans-relativistic plasma. This 

process results in a power-law spectrum extending to energies 

determined by the electron temperature in the hot “corona” (for a 

detailed discussion, see Rybicki & Lightman 1979). The power- 
law index is a function of the plasma temperature, T, and optical 

depth, r. This scenario describes the hard X-ray/soft y-ray 

spectra of bright Seyferts reasonably well (see, e.g., Zdziarski 

et al. 2000). 

There are few physical constraints on the nature of the corona. 
Broadband UV /X-ray spectra require that it not fully cover the 

disk, and suggest it is probably patchy (Haardt et al. 1994). 
X-ray microlensing experiments suggest it is compact; in some 
bright quasars a half-light radius of r c 6 r g has been measured 

(Chartas et al. 2009; Reis & Miller 2013), where we define r c as 

the radius of the corona and r g = GM/c 2 . Eclipses of the X-ray 

source have also placed constraints on the size of the hard X-ray 

emitting region(s): r c Sj 1 70 r„ (e.g., Risaliti 2007; Maiolino 

et al. 2010; Brenneman et al. 2013). 

However, the coronal temperature and optical depth remain 
poorly constrained due to the lack of high-quality X-ray mea- 

surements extending above lOkeV. Complex spectral compo- 
nents, including reflection from the accretion disk as well as 

from distant matter, contribute to the ~ few — 30 keV spectrum, 


and constraints from previous observations suggest typical cut- 
off energies are £ cut Jj 150 keV, requiring spectra extending 
above 50keV for good constraints. The NuSTAR high-energy 
focusing X-ray telescope (Harrison et al. 2013), which covers 
the band from 3-79 keV with unprecedented sensitivity, pro- 
vides the capability to measure both the Compton reflection 
component from neutral material and to constrain the contin- 
uum cutoffs in bright systems. 

The nearby Seyfert galaxy IC 4329A (z = 0.0161, Willmer 
et al. 1991; Galactic /V H = 4.61 x 10 20 cm -2 , Kalberla et al. 
2005; M bh = 1.20 x 10 8 M 0 , de La Calle Perez et al. 2010) is 
a good candidate for such measurements. With a 2-10 keV flux 
range of F 2-10 ~ (0. 1-1.8) x 10~ 10 ergcnT 2 s _1 (Beckmann 
et al. 2006; Verrecchia et al. 2007), IC 4329A is one of the 
brightest Seyferts. In the hard X-ray/soft y- ray band it appears 
similar to an average radio-quiet Seyfert (e.g., Zdziarski et al. 
1996). A Compton reflection component and strong Fe Ka line 
(Piro et al. 1990) are both present. 

IC 4329A has been observed in the hard X-ray band by 
BeppoSAX (Perola et al. 2002), ASCA+RXTE (Done et al. 2000), 
and INTEGRAL (Molina et al. 2013), which placed rough con- 
straints on the high-energy coronal cutoff at £ cut 7s 180keV, 
E cut = 150-390 keV, and E cut = 60-300 keV, respectively. 
Combining INTEGRAL and XMM-Newton data further con- 
strained the cutoff energy to £ cut = 152^ 5 3 1 2 keV (Molina et al. 
2009). The soft X-ray spectrum is absorbed by a combina- 
tion of neutral and partially ionized gas, with a total column 
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of ~3 x 10 21 cm~ 2 (Zdziarski 1994; Madejski et al. 1995; 
Steenbrugge et al. 2005), comparable to the host galaxy’s inter- 
stellar medium (ISM) column density (Wilson & Penston 1979). 
After accounting for the reflection component, the intrinsic spec- 
tral variability is modest (Madejski et al. 2001 ; Miyazawa et al. 
2009; Markowitz 2009). 

We report here on the NuSTAR portion of our simultaneous 
NuSTAR+Suzaku observation of IC 4329A. The combined 
NuSTAR+Suzaku analysis will be described in a forthcoming 
paper (L. W. Brenneman et al., in preparation). In Section 2, 
we report on the NuSTAR observations. Our spectral analysis 
follows in Section 3, with a discussion of the coronal properties 
and their implications in Section 4. 

2. OBSERVATIONS AND DATA REDUCTION 

IC 4329A was observed by NuSTAR quasi-continuously 
from 2012 August 12-16. After eliminating Earth occultations, 
passages through the South Atlantic Anomaly and other periods 
of high background, the NuSTAR observation totaled ~162ks 
of on-source time for focal plane module A (FPMA) and 
159 ks for focal plane module B (FPMB). The total counts after 
background subtraction in the 3-79 keV band were 426, 274 
and 403, 588 for each instrument, yielding signal to noise ratios 
(S/Ns) of 26.8 and 24.6, respectively. 

The NuSTAR data were collected with the FPMA and FPMB 
optical axes placed roughly ~2 arcmin from the nucleus of 
IC 4329 A. We reduced the data using the NuSTAR Data Analysis 
Software (nustardas) and calibration version 1.1.1. 13 We 
filtered the event files and applied the default depth correction 
using the nupipeline task. We used circular extraction regions 
75 arcsec in radius for the source and background, with the 
source region centered on IC 4329A and the background taken 
from the corner of the same detector, as close as possible to the 
source without being contaminated by the point spread function 
wings. Spectra, images and light curves were extracted and 
response files were generated using the nuproducts task. In 
order to minimize systematic effects, we have not combined 
responses or spectra from FPMA and FPMB, but instead fit 
them simultaneously. We allow the absolute normalization for 
both modules to vary, and we find a cross-calibration factor of 
1 .072 ± 0.002 for FPMB relative to FPMA. 

For all the analysis we usedxsPEC version 12.8.1, along with 
other ftools packages within HEASoft 6.14. All errors in the 
text are I o confidence unless otherwise specified, while the final 
parameter values and their uncertainties are quoted in Table 1 at 
90% confidence. 

IC 43 29 A demonstrated a modest, secular flux evolution 
during our observation, increasing by ~ 1 2% over the first 
50 ks of the observation (using clock time), plateauing at 
maximum flux for ~50ks then decreasing by ~34% over the 
remainder of the observation. On average, the flux we measure, 
F 2-10 = 1 -02 x I O' 10 erg cm~ 2 s' 1 , is within its historical range 
of values. Given the lack of short term variability and modest 
flux evolution we use the time-averaged spectrum for all spectral 
analysis (Section 3). 

3. SPECTRAF ANAFYSIS 

The focus of this paper is the characterization of the primary 
high-energy continuum in IC 4329A, and we therefore restrict 
our analysis to 5-79 keV to avoid the effects of low-energy 


13 http://heasarc.gsfc.nasa.gov/docs/nustaiV 


absorption which are poorly constrained by NuSTAR alone. We 
assess the contributions of distant reflection from the outer disk 
and/or torus, although we defer a detailed discussion of the 
distant and inner disk reflection to a forthcoming paper (F. W. 
Brenneman et al., in preparation). 

3.1. Phenomenological Modeling of the NuSTAR Spectra 

We first consider a phenomenological model for the NuSTAR 
spectra to generally investigate the prominence of the high- 
energy cutoff. We begin by evaluating the 5-79 keV data against 
the pexmon model (Nandra et al. 2007), which incorporates 
both primary emission (in the form of a power law with an 
exponential cutoff) and distant reflection. Pexmon is based on the 
pexrav model of Zdziarski et al. (1995), but in addition to the 
reprocessed continuum emission it also includes the fluorescent 
emission lines expected to accompany the Compton reflection 
(Fe Ka, Fe K/3, Ni Ka, and the Fe K« Compton shoulder). 

Fitting the 5-79 keV data using pexmon with no exponential 
cutoff or reflection ( E cui = 10 6 keV and R = 0) yields 
a goodness-of-fit of x 2 / y = 2704/1723(1.57), with clear 
residuals remaining in the Fe K band (Figure 1, top). The 
power law has a slope of T = 1.71 ± 0.01. The spectrum 
has a pronounced convex shape characteristic of both Compton 
reflection and a high-energy cutoff above ~30 keV. 

Allowing the reflection component to fit freely, we get 
R = —0.42 ± 0.02 (negative because of the way reflection 
is parameterized within the model; the absolute value is the 
real reflection fraction) with T = 1.83 ± 0.02, for x 2 /v — 
2010/1722(1.17). Clear residual curvature above ~25keV 
remains (Figure 1, middle). Also freeing the cutoff energy 
of the primary continuum yields £ cut = 149±16keV with 
r = 1.70 ±0.02 and R = -0.34 ± 0.03 with / 2 /v = 
1881/1721 (1.09) (Figure 1, bottom). It is clear that both 
Compton reflection and a high-energy cutoff are required: the 
improvement in fit of Ax 2 /Au = —129/ — 1 (~7%) upon 
addition of a high-energy cutoff is highly significant. 

Freeing the iron abundance (while keeping the abundances of 
other elements fixed to their solar values) improves the fit only 
slightly to y 2 /v = 1871/1720(1.09), with Fe/solar= 1.57 ± 
0.22, though the uncertainties on the other parameters increase 
slightly as a result. Allowing the inclination angle of the reflector 
to fit freely yields no constraints on the parameter and no further 
improvement in fit, so we have elected to keep it fixed at i = 60° . 
We note that when the pexmon component is replaced with the 
more common model of pexrav plus Gaussian emission lines, 
the fit yields similar values of T, R, and E cat . 

Residuals still remain in the Fe K band, suggesting the 
presence of an underlying broad component of the Fe Ka 
line. When this feature is modeled with a Gaussian ( E = 
6.51 ± 0.05 keV, a = 0.36 ± 0.03 keV, EW = 60 ± 15 eV), 
the goodness-of-fit improves to / 2 /v = 1831/1717(1.07). 
Including this component also lowers the iron abundance to 
Fe/solar = 0.87 ± 0.21. We refer to this model hereafter as 
Model 1 . The modest broad iron line detection will be discussed 
at length in F. W. Brenneman et al. (in preparation), and will 
not be further addressed in this work. 

3.2. Toward a More Physical Model 

The high S/N of the data enable us to consider more 
physically motivated models to describe the continuum emis- 
sion. We employ the compPS model of Poutanen & Svensson 
(1996), which produces the continuum through inverse Compton 
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Figure 1. Plots of the data-to-model ratio of the NuSTAR/FFMA (black) and FPMB (red) spectra vs. a pexmon model without an exponential cutoff or reflection (top 
left), including reflection (top right), including an exponential cutoff in addition to reflection (bottom left), and including also a broad Gaussian Fe Kc*f line (bottom 
right). The horizontal green line in each plot represents a perfect fit. 

(A color version of this figure is available in the online journal.) 


Table 1 

Best-fit Parameters, their Values and Errors (to 90% Confidence) for Models 1-3 


Component 

Parameter (units) 

Model 1 

Model 2 

Model 3 

TBabs 

IVh (x 10 20 cm -2 ) 

4.61(/) 

4.61(/) 

4.61(/) 

pexmon 

r 

1-70- 0 o°o 4 4 

1.70(/) 

1-70 (/) 


Ecnt (keV) 

178 + J 4 o 

90^8* 

96*f 4 * 


Kpex (x 10 -2 photons cm -2 s -1 ) 

2 90 +0 - 20 
z ‘ yu -0.18 

2.90(f) 

2.90 (/) 


R 

n 9Q+0.06 
-u.zy_ 006 

—0 33 +0 - 04 
u - :o -0.04 

-0-32To4 


Fe/solar 

O.87 + _° 0 4 3 2 

0 63 +0 20 
u -° -0.20 

°' 68 -0 2 19 

compTT 

kT e ( keV) 


33*6* 

37*6* 


r 


3-41 40 o 5 3 8 8 

i oc+0.20 
-0.10 


Kcom (x 10 -3 photons cm -2 s _1 ) 


9 Q 1+1 .88 

V.ui_ 29 5 

o i 7+2.37 
81 —2.96 

zgauss 

Eb (keV) 

6-51 + -°o°05 

f\ 49+0.07 
oo6 

A8+0.08 

b-45_oo6 


a b ( keV) 

0-36- 0 0 12 

0 36 +0 09 
U "- )O -0.08 

Q 77+O.IO 
yjmD -0.08 


K b (x 10~ 5 photons cm~ 2 s _1 ) 

7-84-2 7 6 2 l 

9 09+2-59 
y - uy -2.50 

q i i +2.64 
y i l -2.54 


EW b (eV) 

60^0 

70+20 
/u — 19 

7Q+20 

/u -20 

Final fit 


1831/1717(1.07) 

1836/1718(1.07) 

1837/1718 (1.07) 


Notes. Parameters marked with an (f) are held fixed in the fit, while those marked with an (*) are tied to another parameter. 
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scattering of thermal disk photons off of relativistic electrons sit- 
uated above the disk in a choice of geometries. The disk photons 
and coronal plasma are each parameterized by a single temper- 
ature. It then reflects the continuum X-rays off of a slab of gas 
(the disk or torus) to produce reprocessed continuum emission 
self-consistently by internally linking with the pexrav model. 
Reprocessed line emission is not included. 

We fix the energy of the thermal disk photons to kT^ = 
30 eV, appropriate for a black hole of M B h ~ 10 s M 0 (Frank 
et al. 2002). We initially fix the geometry of the hot electrons to 
be spherical with a Maxwellian distribution, the former choice 
being based on the coronal compactness measurements cited in 
Section 1 (e.g., Chartas et al. 2009), although we also compare 
the results with a slab geometry. We fit for the electron energy 
( kT e ) and coronal optical depth (r), as well as for the reflection 
fraction (R) of the reprocessing gas, and model normalization 
( K ). We add in separate Gaussian components to represent the 
Fe K a (narrow and broad) and K/J line emission resulting from 
reflection in order to maintain consistency with the pexmon 
model employed in Section 3.1 above. 

This approach yields approximately the same global 
goodness-of-fit as Model 1: / 2 /v = 1849/1715(1.07), with 
parameter values of kT e = 48 ± 22keV (assuming that E cut = 
3kT e , as per Petrucci et al. 2001, this is equivalent to a cut- 
off energy of approximately 144 ± 66keV, consistent with that 
measured in Model 1), r = 2.70± 1.31 and R = — 0.39±0.06, 
assuming a distant, neutral reflector inclined at 60° to the line of 
sight. We note that the iron abundance is not constrained by the 
model. The narrow Gaussian components representing narrow 
Fe Ktf and K/i have equivalent widths of FW Ka = 41 ± 8eV 
and EWj^i = 15 ± 12 eV, respectively. The addition of these 
components improves the global goodness-of-fit by A/ 2 = 483 
and A/ 2 =19, respectively, each for three additional degrees 
of freedom. 

We also note that consistent results are achieved (within er- 
rors) when we fit for the Compton-y parameter rather than 
the optical depth, as per the approach taken in Petrucci et al. 
(2013): kT e = 43 ± 41 keV and y = 1.16 ± 0.05, resulting 
in r = 3.45 ± 0.95. In addition to having larger parameter 
uncertainties, however, this approach also results in an uncon- 
strained reflection fraction. Consistent results are obtained when 
we fix the reflection fraction at R = —0.39 (as described in the 
previous paragraph) and Fe/solar = 1: k 7/ = 39± 9keV, 
y = 1.06 ± 0.01 and r = 3.47 ± 0.25. However, in light of 
the importance of probing the reflection as a free parameter, we 
have elected to fit for the optical depth explicitly in the model. 

We consider a slab geometry for the corona as well, again 
with a Maxwellian electron distribution. This results in no 
significant improvement to the fit: x 2 / y = 1836/1713(1.07), 
with kT e = 48 ± 13keV (equivalent to a cutoff energy of 
£ CU [ = 144 ± 39keV, also consistent with that measured in 
Model 1), t = 1 .50 ± 0.39 and R = -0.38 ± 0.07. This model 
is also insensitive to the Fe abundance. The equivalent widths 
of the Gaussian lines are the same as for the spherical geometry, 
within errors. 

For both geometries the electron plasma temperature is 
constant within errors, while the optical depth is pushing 
the upper limit of its sensible parameter space. We therefore 
check the compPS results using the compTT model (Titarchuk 
1994), which assumes a simpler thermal electron distribution 
and does not self-consistently include the Compton reflection 
continuum. We add in reflection from distant matter using 
the pexmon model. We assume an incident power-law photon 
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Figure 2. Upper panel: EFe vs. energy plot for Model 1 vs. the FPMA (black 
points) and FPMB (red points) data. The pexmon model incorporating the cutoff 
power-law continuum and distant reflection is in dashed red, the direct power- 
law emission is in dashed blue, and the broad Gaussian in dashed green. The 
summed models are in solid black (FPMA) and solid red (FPMB). Lower panel: 
data-to-model ratio for the FPMA (black) and FPMB (red). The horizontal green 
line shows a ratio of unity for reference. 

(A color version of this figure is available in the online journal.) 

index of T = 1.70 (fixed), as found in Model 1, and tie 
the normalization of the reflected emission to that of the 
Comptonized component, such that the contribution of the 
reflector relative to the Comptonization is entirely determined by 
R and the iron abundance. We checked the models to ensure a 
tight match between the shapes of a power law of this index 
and the shape of the compTT component. We also fix the 
cutoff energy of the incident power law at E cut = 3 kT e . The 
spherical geometry compTT component model and the same 
model with a slab geometry are referred to hereafter as Models 
2-3, respectively. 

We find values for the coronal temperature consistent within 
errors with those of compPS for both the spherical and slab 
geometries: kT e = 33 ± llkeV and k 7/ = 37 ± 16keV, 
respectively. The coronal optical depth in the slab geometry is 
slightly smaller with compTT versus compPS, though consistent 
within errors: r = 1.25 ± 0.50. For the spherical geometry, 
r = 3.41 ± 0.90, also consistent with its compPS analog within 
errors. The reflection fraction measured with this model (defined 
in the same way for both geometries) is comparable to that 
determined by compPS: R = —0.33 ± 0.03 for the sphere 
and R = —0.32 ± 0.04 for the slab. The iron abundance is 
comparable that that found in Model 1: Fe/solar = 0.63 ± 0.20 
(sphere) and Fe/solar = 0.68 ± 0.20 (slab). 

Table 1 provides the best-fit parameters and their Markov 
Chain Monte Carlo (MCMC) derived 90% confidence errors 
for all three models. Figure 2 plots the contributions of in- 
dividual model components to the overall fit. We show only 
the components for Model 1, since Models 2-3 look vir- 
tually identical, except using a Comptonization component 
in lieu of a power law. The total absorbed 5-79 keV flux 
and luminosity are F5-79 = 3.04 x 10~ 10 erg cm" 2 s -1 and 
L5-79 = 1.77 x 10 44 ergs~' , respectively. 

The MCMC analysis employed to determine the formal 
parameter distribution used the Metropolis-Hastings algorithm 
(e.g., Kashyap & Drake 1998 and references therein) following 
the basic procedures outlined in, e.g., Reynolds et al. (2012). 
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Figure 3. Output from the MCMC analysis of Model 1: contours show the 
67% (black), 90% (red) and 99% (green) probability densities for R vs. E cnt . 

(A color version of this figure is available in the online journal.) 

Using this basic procedure we generated probability density 
contours for the most interesting pairs of parameters for each 
model, shown in Figures 3 and 4. Both sets of contours are 
closed, implying that kT e and r are independently constrained 
to 99% confidence in both the spherical and slab geometries. 
Nonetheless, some degeneracy between the two parameters 
still remains, as evidenced by the linear correlation seen in 
each plot, due to an inherent modeling degeneracy between the 
optical depth and temperature of the electron plasma in each 
geometry. Figure 4 depicts the modest range over which these 
parameters are degenerate, which is ±18% of the parameter 
space in temperature and + 1 ,™ in optical depth for the sphere, 
versus in temperature and in optical depth for the slab. 

4. DISCUSSION 

The time-averaged 5-79 keV NuSTAR spectrum is well- 
described by a largely phenomenological pexmon model 
(Model 1) that includes, to 90% confidence, a continuum 
power law with F = 1 .70 ± 0.03, and a high-energy cutoff 
of £ cut = 178^40 keV, as well as reprocessed emission from 
a distant reflector ( R = —0.29^°^, Fe/solar = 0.87*0 3 ^). 
Prior measurements provided constraints on the cutoff energy 
of the power law at E cut ^ 180keV (Perola et al. 2002), 
E c at = 1 50-390 keV (Done et al. 2000), and L cul = 60-300 keV 
(Molina et al. 2013), though Molina et al. (2009) did achieve 
a more precise constraint of E cul = 1 52 + 32 keV using XMM 
in tandem with INTEGRAL. The two datasets were not taken 
simultaneously, however, and the S /N achieved by the NuSTAR 
data is superior to that of INTEGRAL. We therefore consider 
our new measurements — which agree with all of the previously 
mentioned values, within errors — to be more robust. 

The high S /N, simultaneously obtained broad band NuSTAR 
spectrum enables us to apply physical models for the underlying 
coronal continuum emission that go beyond phenomenological 
descriptions. The models parameterize the temperature and 
optical depth of the electron plasma for two coronal geometries; 
a sphere and a slab. Both geometries fit the data equally well, 
though we note that the sphere model provides slightly tighter 
constraints on several of its parameters. Both models also 
produce consistent values for the electron temperature within 
errors. Flowever, they result in slightly different values for the 



kT e (keV) 



Figure 4. Output from the MCMC analysis of Model 2 (sphere; upper panel) 
and Model 3 (slab; lower panel): contours show the 67% (green), 90% (red) and 
99% (black) probability densities for r vs. kT e . 

(A color version of this figure is available in the online journal.) 

optical depth: kT e — 37+ 7 6 keV with r = 1.25^° ( j 2 I 0 0 for the slab 
geometry, compared with kT e — 33+ 6 6 keV with r = 3.4 Uq 5 3 8 8 
for the spherical geometry (both at 90% confidence). 

This discrepancy in optical depth is due primarily to the 
way that the value is calculated for a given geometry within 
the Comptonization models we employ (e.g., Titarchuk 1994; 
Poutanen & Svensson 1996): the optical depth for a slab 
geometry is taken vertically, whereas that for a sphere is taken 
radially and thus incorporates an extra factor of 1 /cos(60) = 2. 
Taking this extra factor into account, the optical depth of the 
spherical case can be translated into the slab geometry for ease of 
comparison: r = 1.7 1^ 2 ® for the sphere versus r = 1.25+o 2 0 
for the slab. Though these values do not formally agree within 
their 90% confidence errors, they are compatible at the 2a level. 

The derived electron temperatures for the sphere and slab 
coronal geometries are low compared to E c ut /2, but are not far 
off from £ cut /3, which is consistent with the corona having 
significant optical depth (i.e., r > 1), modulo uncertainties in 
geometry, seed photons, outflows, anisotropy, etc. which we 
are not able to probe even with our high S /N data. Due to an 
inherent modeling degeneracy between the optical depth and 
temperature of the electron plasma in each geometry, there is a 
small, linearly correlated range of values for these parameters 
which demonstrate approximately equal statistical fit quality. 
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as can be seen in Figures 3 and 4. Nonetheless, we constrain 
both parameters precisely and accurately with the best data ever 
achieved over this energy band. The data quality and goodness- 
of-fit of the Comptonization models gives us confidence in the 
temperatures and optical depths we have measured. We note, 
however, that without high S/N data at energies ^lOOkeV we 
are unable to discriminate between a thermal and non-thermal 
population of coronal electrons (e.g., with a model such as 
eqpair, Coppi 1999). A significant non-thermal contribution 
could change the temperatures and optical depths that we 
measure. 

With our robust determination of the continuum shape over 
the broad energy range, we estimate that the power dissipated in 
the corona, in the form of the power-law continuum, is ~87% 
of the total luminosity of the entire system from 5-79 keV (the 
power law has a luminosity of Ls _ 79 = 1.53 x 10 44 ergs' 1 ). 
This represents ~ 10% of the bolometric luminosity of the source 
(Lboi = 1.60 x lO^ergs' 1 ; de La Calle Perez et al. 2010). 
Given that IC 4329A is a radio-quiet active galactic nucleus 
(total 10 MHz-100 GHz L r = 2.2 x I C 39 erg s ' 1 ; Wilson & 
Ulvestad 1982), we do not expect a significant portion of the 
X-rays to come from a jet component, so we may infer that 
the remaining ~13% of the 5-79 keV emission represents the 
contribution from reflection. 

Our spectral fitting results are broadly consistent with the sig- 
natures expected from dynamic, outflowing coronae as defined 
by Beloborodov (1999) and Malzac et al. (2001): a hard spec- 
tral index and relatively weak reflection. If the corona is really 
powered by compact magnetic flares that are dominated by e ± 
pairs, the resulting plasma is subject to radiation pressure from 
the photons moving outward from the disk. The bulk velocity 
at which the plasma should be outflowing from the disk, given 
our Model 1 photon index of T = 1.70, is i> = 0.22c, accord- 
ing to Beloborodov (1999). Our measured reflection fraction 
of R — 0.3 is also consistent with that predicted by Malzac 
et al. (2001) for a Comptonizing plasma with these parameters. 
We note that the beaming of the coronal emission inherent in 
such dynamic models not only hardens the spectrum, but also 
implies that the intrinsic break in the spectrum occurs at lower 
energies than what is observed: i.e., £ cut ~ 100 keV rather than 
£cut = 178 4 4 q keV, as measured with Model 1 . A lower intrin- 
sic rollover is consistent with our measurements of the electron 
temperature from Comptonization models ( kT ~ 35keV), if 
we assume that E cut is between 2-3 kT e . 
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